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r- H ■ Abstract. Deep VLT optical spectroscopy, HST+ACS (GOODS) imaging and VLA observations are used to 

unveil the nature of a complete sample of 47 EROs with R — K s > 5 and K s < 20. The spectroscopic redshift 
cm : completeness is 62%. Morphological classification was derived for each ERO through visual inspection and surface 

brightness profile fitting. Three main ERO morphological types are found: E/S0 galaxies (~ 30 — 37%), spiral-like 
(~ 24 — 46%) and irregular systems (~ 17 — 39%). The only ERO detected in the radio is likely to host an obscured 
AGN. The average radio luminosity of the star- forming EROs undetected in the radio implies star formation rates 
of the order of ^33 M0yr -1 . The colors, redshifts and masses of the E/S0 galaxy subsample imply a minimum 
formation redshift Zf ~ 2. With this Zf there is enough time to have old and massive stellar spheroids already 
^ ■ assembled at z ~ 1. We verify that the R — K s vs. J — K s color diagram is efficient in segregating old and 

^ dusty- star- forming EROs. 
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1. Introduction challenging even with 8-10m-class telescopes. However, as 

different morphologies are expected in case of spheroids or 
Extremely Red Objects (EROs) are important probes of star _ forming syste ms, HST imaging can be used to corn- 
galaxy evolution. On one hand, they allow to constrain the plement spectroscopy for the faint ERO population (e.g. 
- formation epoch, the evolutionary pattern andtfie clus- Moriondo et aL 2000; Yan & Thompson 2003). The main 
tering of the oldest early-type galaxies at z ~ 1 (e.g. open questions about EROs include: the relative fraction 
Daddi et al. 2000a,b; McCarthy et al. 2001; Cimatti et f different ERO types, the nature and star formation 
al. 2002a; Smith et al. 2002; Firth et al. 2002; Moustakas rates of dusty EROs, the number density of the oldest 
& Somerville 2002; Roche et al. 2002, 2003; Saracco et n igh-z spheroids selected as EROs, the fraction of ob- 
al. 2003). On the other hand, they allow to select dusty S cured AGN hosted by EROs. In this paper we present the 
star-forming galaxies and obscured AGN in a way com- main results of very deep VLT spectroscopy, HST imaging 
plementary to submillimeter surveys (Cimatti et al. 1998; anc j VLA observations aimed at investigating the nature 
Afonsoetal. 2001; Mohan et al. 2002; Cimatti et al. 2002a; Q f a complete sample of 47 EROs with K 8 < 20. The 
Franceschini et al. 2002; Brusa et al. 2002). Although properties of the sample will be presented in a companion 
recent VLT spectroscopy shed new light on moderately p ape r (Cassata et al. 2003, in preparation). H = 70 km 
bright EROs with K s - 19 (Cimatti et al. 2002a), the s" 1 Mpc" 1 , ft m = 0.3 and ^ A = 0.7 are adopted, 
characteristics of fainter EROs is still unclear mainly be- 
cause of their weakness which makes spectroscopy very 

2. The sample and the observations 

Send offprint requests to: Andrea Cimatti, e-mail: 

cimatti@arcetri.astro.it The complete sample (already described by Cimatti et al. 

* Based on observations made at the European Southern 2002a) is made by the 47 EROs with R - K s > 5 and 

Observatory, Paranal, Chile (Program 70.A-0140(B)). K s < 20 present in the 32.2 arcmin 2 area of the Chandra 
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Deep Field South (CDFS; Giacconi et al. 2001) covered 
by the K20 survey (Cimatti et al. 2002b). One object was 
removed from the original sample because it was found to 
be a blend of two objects both with K s > 20. 

New VLT-UT4+FORS2 multi-object optical spec- 
troscopy was made in November 2002 with the red- 
sensitive CCD array. The total integration time was 7.8 
hours during 0.7"-1.0" seeing conditions, with 1.0" wide 
slits, grism 3001 (R~660) and dithering of the targets 
along the slits. The mask included 13 EROs of the 
K20/CDFS sample that were previously unidentified or 
not observed. The spectra were reduced and analyzed 
with the methods described Cimatti et al. (2002a, 2002b). 
VLT spectroscopy was complemented with the HST+ACS 
imaging in the BViz bands taken in the GOODS/HST 
Treasury Program (Giavalisco et al. 2003). 

Deep radio data were obtained with the VLA at 1.4 
GHz in 1999 October with a bandwidth of 25 MHz in 
each of two circular polarizations, and at 4.8 GHz in 1999 
Feb with a bandwidth of 100 MHz in each of two circu- 
lar polarizations, for 12-hours total at both frequencies 
(Kellermann et al. 2003, in preparation). The rms sen- 
sitivities were 8.6 /iJy and 9.6 /iJy at 1.4 and 4.8 GHz, 
respectively, each with a full-width half-power resolution 
of 3.5". The reductions were made using AIPS and each 
field was self-calibrated. The position and flux density of 
each source with a peak flux density greater than 3.5-rms 
was determined from a Gaussian model fit of the source. 

3. VLT spectroscopy and HST imaging results 

The 13 ERO spectra provided 9 new redshift mea- 
surements: 7 star- forming galaxies at z ~ 1.2 — 1.6 
with [OIIJA3727 emission, and 2 early- type objects at 
z ~ 1.5. When added to the previous identifications of 
Cimatti et al. (2002a), the spectroscopic redshift com- 
pleteness for the EROs in the K20/CDFS area increases 
to 81% (21/26) (K s < 19.2) and 62% (29/47) (K s < 
20.0). Photometric redshifts for the 18 spectroscopically 
unidentified EROs were derived using the ESO/ GOODS 
VLT+FORS1 BVRIz and the new VLT+ISAAC JHK S 
public imaging and the Hyper z software (Bolzonella et al. 
2000). The agreement between spectroscopic and photo- 
metric redshifts is good and does not depend on the ERO 
class, with Az — (z spec — %>hot)=0.001 and <r(Az)=0.16 
(29 EROs). Figure 1 shows the redshift distributions di- 
vided into morphological classes (see Section 4). 

The morphological type of EROs was first estimated 
by visual inspection of the ACS images and isophotal 
maps in the four available bands (BViz). One ERO (not 
observed spectroscopically) is unresolved and has colors 
consistent with being a M star. Since in some cases the 
spatial resolution and/or the S/N ratio turned out to 
be limited, we just assigned each object to one of the 
following broad classes: E/S0 galaxies (class 1), spiral- 
like objects (class 2) and irregular /interacting systems 
(class 3). Ambiguous cases were flagged by introducing 
the intermediate classes 1.5 and 2.5. Thus, we used the 
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Fig. 1. The redshift distributions. The shaded histograms in- 
dicate the photometric redshifts. 



GASPHOT software (Pignatelli and Fasano 1999) to carry 
out a quantitative analysis on the z-band ACS images 
(filter F850LP) by fitting the surface brightness profile 
of each ERO with the Sersic law (SB oc r 1//n ), convolved 
with the average PSF extracted from the stars in the field. 
Being the bulk of the sample at z ~ 1, the filter F850LP 
is close to the rest-frame 5-band, thus minimizing the ef- 
fects of dust extinction and morphological K-correction 
for the available set of ACS filters. EROs belonging to the 
class 1.5 and having n >2.5 or <1.5 were assigned to the 
class 1 or 2, respectively, while the cases with 1.5< n <2.5 
were assigned to intermediate class 1.5. Instead, the Sersic 
index showed to be ineffective in disentangling classes 2 
and 3, making advisable to leave unchanged most of the 
dubious cases with flag 2.5 (10 objects). Indeed, for six 
of these objects the fitting algorithm failed to converge, 
due to the extremely irregular and noisy image structure. 
These objects are not reported in the lower panel of Fig. 2, 
where the final classification is plotted versus the Sersic in- 
dex. We note that the larger number of class-2.5 vs. class- 
1.5 objects and the higher redshifts of "irregular" EROs 
(Fig.l) may be due to morphological K-correction effects 
(e.g. Marcum et al. 2001; Windhorst et al. 2002). By tak- 
ing into account the previous uncertainties, the ranges of 
morphological fractions are 30-37% (class 1 - class 1+1.5), 
24-46% (class 2 - class 2+2.5) and 17-39% (class 3 - class 
2.5+3). The fractions are in reasonable agreement with the 
recent results by Moustakas et al. (2003) and Gilbank et 
al. (2003), and do not significantly change as a function 
of the limiting K s magnitude. Fig. 2 (top panel) and Fig. 3 
show that a rather good agreement is present between the 
morphological and spectroscopic classes. 
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Fig. 2. The Sersic index n vs. the spectral (top) and visual 
morphological class (bottom). Filled symbols: reliable fit of n; 
open symbols: uncertain fit of n. Spectral classes: -1 = uniden- 
tified or unobserved, 1 = passive early- type, 1.5 = early- type 
+ weak [OII]A3727 emission, 2 = star- forming. Visual morpho- 
logical classes: 1 = E/SO, 2 = spiral- like, 3 = irregular. 



4. Dusty star-forming EROs and obscured AGN 

Only one object was detected at > 5<r level (S > 50 
/iJy) in the 1.4 GHz radio source catalog of the CDFS 
(Kellermann et al. 2003, in preparation). Our radio de- 
tection rate (1I0J3/46 at la poissonian confidence level) 
is compatible with the preliminary results of the Phoenix 
Deep Survey (Hopkins et al. 2003) where ~6% of a sam- 
ple of -425 EROs (K s < 20, R - K s > 5) are detected 
with fluxes > 60 /iJy (5a) (J. Afonso, private commu- 
nication), with the results of Smail et al. (2002) (3/50 
EROs with K 8 < 20, R - K s > 5.3 and S > 50 /iJy at 1.4 
GHz) and Roche et al. (2002) (1/31 EROs with K s < 19.5, 
R - K s > 5 and S > 50 /iJy at 1.4 GHz). 

The ERO detected in the radio is an emission line 
object with [OIIJA3727 at z = 1.610 and a spiral-like 
morphology. The flux densities at 1.4 GHz and 4.8 GHz 
are 106.8±8.6 fiJy and 155.1±9.6 /iJy respectively, im- 
plying a flat /inverted spectrum with a ra dio = +0.3 
(S(u) oc v a ). The object is also present in the X-ray cat- 
alog of Alexander et al. (2003) (RA=03:32:10.80, Dec= 
-27:46:27.6, J2000) with fluxes F(0.5-2 keV)=8.67x lO" 17 
cgs and F(0.5-8 keV)=3.08x 10 -16 cgs, corresponding to 
luminosities of 1.4xl0 42 erg s _1 and 5xl0 42 erg s _1 re- 
spectively. Its X-ray luminosity, owio and location in typ- 
ical segregation diagrams (e.g. a ra di vs. radio flux, optical 
flux vs. 0.5-8 keV flux; Bauer et al. 2002) are more sugges- 
tive of a Type 2 AGN than a starburst galaxy. The lack 
of broad MgIIA2800 emission implies obscuration of the 
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Fig. 3. Examples of EROs VLT spectra and HST morpholo- 
gies (ACS filter F850LP, box size = 4x4 arcsec). From top: 
E/S0 (z = 1.096), spiral- like (z = 1.222), irregular (z = 1.221). 



active nucleus. From the X-ray hardness ratio we derive 
log(NH)< 22.7 cm -2 (assuming a photon index T = 1.7). 

Although EROs are often found as faint counterparts 
of X-ray selected sources (e.g. Yan et al. 2003), the fraction 
of if -selected EROs with X-ray emission is rather small 
(e.g. Alexander et al. 2002; 14^% at K <20.1; Roche et 
al. 2003; -8% at K s < 20). Taking into account the only 
other X-ray source detected in our ERO sample (Brusa et 
al. 2002), the fraction of X-ray emitting EROs amounts 
to 2/46 (1.5-10% at la poissonian confidence level), con- 
sistently with the above mentioned results. 

Although only one ERO has been detected signifi- 
cantly in the radio, a mean flux density of the undetected 
sources was estimated by comparing the flux density dis- 
tribution of the non-detected ERO sources with that ex- 
pected from noise and confusion. The radio intensity dis- 
tribution at 1.4 GHz for the 28 non-detected spiral-like + 
irregular ERO sources has a rms of 9.3 /iJy and an aver- 
age value of 3.5 /iJy compared with an rms of 8.9 /iJy and 
average value of —0.8 /iJy for the control sample of 270 
random locations in the field. From this excess we infer 
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that the average flux density of ERO sources at 1.4 GHz 
is 3.5 ± 1.7 /iJy. At 5 GHz we estimate an average flux 
density of < 4 jaJy. For the 17 E/SO galaxies we obtained 
limit of < 3 /iJy at both 1.4 and 5 GHz. 

Adopting a radio spectral index a = —0.6, the av- 
erage rest-frame 1.4 GHz luminosity is 2.7 x 10 22 W 
Hz -1 for the average redshift of the 28 undetected star- 
forming EROs ((z) = 1.337). Following Condon (1992) 
and Haarsma et al. (2000) (with Salpeter IMF and 0.1-100 
M mass range), we derive SFR~ 33 M yr _1 . The esti- 
mated SFRs are consistent with those obtained by Brusa 
et al. (2002) using the X-ray emission (5-44 M yr _1 ). 
Following Kennicutt (1998), the expected average far-IR 
luminosity is ~ 1 — 2 x 10 11 L corresponding to 850/im 
observed fluxes below 1 mJy, thus explaining the low de- 
tection rate of K s < 20 EROs with submm/mm contin- 
uum observations (e.g. Mohan et al. 2002). In compari- 
son, the [OIIJA3727 luminosity of all the spectroscopically 
identified star-forming EROs in the CDFS provides an av- 
erage SFR=2.8 M yr _1 . If this difference with the radio 
is solely due to dust extinction, the required average red- 
dening is E(B — V) ~ 0.5. This amount of reddening is 
significantly larger than in Lyman-break galaxies (median 
E(B - V) ~ 0.17, range E(B - V) ~ 0.1 - 0.3; Giavalisco 

2002 and references therein). The average SFR(radio) im- 
plies a cosmic SFR density of EROs at z ~ 1.3 of >0.008 
M yr _1 Mpc -3 (adopting the comoving volume between 
^min = 0.726 and z max = 2.185), corresponding to about 
>12% of the total budget at that redshift (e.g. Somerville 
et al. 2001). This is strictly a lower limit because our sam- 
ple is limited to K s < 20 and due to the conservative 
volume adopted between z m i n and z max . 

5. The oldest and massive spheroids at z ~ 1 

The reddest spheroidals selected as EROs represent the 
oldest envelope of the E/S0 population at z ~ 1 and 
are crucial to constrain the earliest epoch of massive sys- 
tem formation. Taking into account the non-poissonian 
fluctuations due to clustering (Daddi et al. 2000a; Roche 
et al. 2002), the observed surface density of E/S0 galax- 
ies (morphological class 1-1.5) is 0.50^o'i6 arcmin -2 (la) 

(0.50to.26 at 2cr )- The Pure Luminosity Evolution (PLE) 
model of Daddi et al. (2000b) with a star formation e- 
folding time r = 0.3 Gyr and a formation redshift zi ~ 
2.5 predicts surface densities of ~2 arcmin -2 and ~0.9 
arcmin -2 for the 2MASS local early-type (E+S0+ some 
contamination by Sa galaxies; Kochanek et al. 2001) and 
the local elliptical (E only; Marzke et al. 1994) luminos- 
ity functions respectively. Although the large uncertainties 
on the surface density prevent a stringent comparison with 
the PLE model, we note that the observed density of E/S0 
galaxies at 0.8 ~ z ~ 1.6 is still formally consistent with 
the above PLE predictions within 1-2<j. The question of 
high-z early-type galaxies is still controversial (e.g. Treu 

2003 and references therein), and large area surveys are 
required to overcome the effects of the strong clustering. 




redshift 

Fig. 4. R - K s vs. z for the EROs with E/S0 morphology. 
Continuous lines indicate the predicted colors for different for- 
mation redshifts (indicated by the numbers) and r = 0.3 Gyr. 

Figure 4 shows a significant spread in the formation red- 
shifts of the identified E/S0 galaxies (2.2 ~ Zf ~ 4.0 for 
r = 0.3 Gyr). If r = 0.1 Gyr, the range is 2 ~ Zf ~ 3. 
From the lowest acceptable formation redshift (zf ~ 2.0) 
to z ~ 1 there is enough time to reach an age of ~2.5 Gyr, 
compatible with that estimated from the average spec- 
trum and from the 4000 A continuum break amplitude of 
old EROs at z ~ 1 (Cimatti et al. 2002a), and within the 
range of the ages of z > 1 spheroidal galaxies (~1.5- 3.5 
Gyr) estimated with different spectral synthesis models 
(e.g. Nolan et al. 2001 and references therein). The rest- 
frame near-IR luminosities of old EROs in our sample are 
in the range of -24.5 and -27 (Pozzetti et al. 2003), and 
their stellar masses are ~ 10 11 -10 12 M (Fontana et al. 
2003, in preparation). 

The morphology, spectral properties, luminosities, in- 
ferred ages, formation redshifts and stellar masses im- 
ply the existence of a substantial population of old, pas- 
sively evolving and fully assembled massive spheroids at 
z ~ 1 — 1.5 which requires that major episodes of mas- 
sive galaxy formation occurred at z ~ 2. Daddi et al. 
(2003) uncovered a population of massive and clustered 
star-forming galaxies at z ~ 2 which may represent the 
long sought-for progenitors of massive spheroids. 

6. Testing the photometric segregation of EROs 

We used our sample to test the segregation of morpholog- 
ically classified EROs in the R — K s - J — K s diagram 
proposed by Pozzetti & Mannucci (2000) (Fig. 5). The 
colors were derived from the new public ESO GOODS 
imaging taken with the VLT FORS1 (R) and ISAAC 
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Fig. 5. The segregation color - color diagram. The E(B-V) 
vector is relative to z — 1.1 and the Calzetti extinction curve. 



fine the boundary line between "old" and "dusty" EROs 
is J — K S =0.34(R — if s ) +0.09. This relation was derived 
specifically for the set of RJK S filters of our photometry, 
and is valid for R-K s >5.3 (Pozzetti & Mannucci 2000). 

Two main results are evident from Fig. 5: (1) the E/S0 
galaxies with secure morphological classification (class=l) 
and most of star-forming EROs lie in their expected re- 
gions, (2) a few star- forming EROs (classes=2-3), how- 
ever, are found in the region of E/S0 galaxies. If we con- 
sider only the EROs with reliable morphological classifica- 
tion, the contamination of the E/S0 region is 4/15 (27%), 
whereas 4/18 (22%) star- forming EROs do not fall in their 
expected locus. The contamination of the E/S0 region may 
be due to star-forming galaxies with a substantial under- 
lying old stellar population which dominates the colors. 
Our results indicate that the color - color segregation is 
efficient and that the typical contamination amounts to 
-20-30% (see also Pierini et al. 2003). 
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